Effects of Andreev reflection on the conductance of quantum-chaotic dots 
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I. INTRODUCTION 

The physics of normal-metal-superconducting (NS) 
hybrid structures has been of interest for several decades. 
In the sixties, the scattering of a quasi-particle at a NS 
interface was first described by Andreev i Thus, it is now 
widely known that an electron approaching to a super- 
conductor from the normal-metal side, with an energy 
above the Fermi level, can be reflected at the interface 
as a hole with an energy below the Fermi level and same 
momentum, but in opposite direction. A Cooper pair 
is thus formed in the superconductor. This scattering 
process is known as Andreev reflection. More recently, 
the interest in this subject was renewed in the field of 
mesoscopic physics, where the electron wave function re- 
mains coherent and therefore quantum interference ef- 
fects play an important role. Several effects of phase- 
coherence Andreev scattering are under current investi- 
gation. In addition, the effects of disorder and ballistic- 
chaotic scattering in hybrid structures have been widely 
studied. For example quantities such as the density of 
states and the Ehrcnfcst time as well as phenomena such 
as coherent backscattering, excitation and gap fluctu- 
ations, superconducting proximity effect, and resonant 
tunneling have been investigated in mesoscopic hybrid 
structures^££i&^±±^^i^ Also, the effects of 
Andreev scattering on the conductance through quan- 
tum dots whose classical dynamics is intcgrable or mixed 
(regular and chaotic) have been investigated. 16 ' 17 

On the other hand, random matrix theory (RMT) 
has been useful to study the statistical properties of 
many different phenomena in normal metal o 18 ' 19 ' 20 ' 21 ' 22 
as well as in superconductors. In particular, the effects 
of NS interfaces on the conductance of disordered and 
chaotic billiards-quantum dots-have received consider- 



able attention^iii^Si 2 !^^^^ Recently, some sta- 
tistical properties of the conductance fluctuations in 
mesoscopic NS heterostructures have been experimen- 
tally studied in Ref. [2j| where the explanation of the ex- 
perimental results was based on RMT of quantum trans- 
port through chaotic quantum dots; it was also pointed 
out the importance of considering the potential barriers 
formed at the NS interface. 

Using a scattering approach to the problem of quantum 
transport and RMT, here we study the statistical prop- 
erties of the conductance through a quantum-chaotic dot 
with a superconducting lead attached to it. In particular, 
we are interested in the distribution of the conductance. 
The configurations of one and two normal leads addition- 
ally attached to the dot are considered. We take into ac- 
count the presence of a Schottky barrier-which is formed 
in real NS junctions- of transparency T in the supercon- 
ducting lead. Thus, for these two setups we calculate the 
complete distribution of conductances for different val- 
ues of r. As we will see, the conductance distribution is 
affected by the presence of the superconducting region as 
well as the scattering at the Schottky barrier. 

This paper is organized as follows. In the next Section 
we briefly review general ideas about electron-hole scat- 
tering at NS junctions. We also introduce the scattering 
matrices associated to the different elements of our ex- 
perimental setup. In Sec. Ill, we study the two-terminal 
configuration. Firstly, we calculate the total scattering 
matrix and the conductance. Secondly, we present our 
results for the distribution of the conductance for the 
cases of time- reversal symmetry ((3 = 1) and broken time- 
reversal symmetry (f3 = 2). Simple analytical expressions 
are provided for the conductance distribution. Sec. IV 
is devoted to the three-terminal configuration. We cal- 
culate the necessary scattering-matrix elements in order 
to calculate the conductance and present the results for 
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its distribution for both symmetries ((3—1 and (3 = 2), 
at different values of the transparency of the Schottky 
barrier. A summary and conclusions of our work are 
presented in Sec. V. 



II. MODEL AND FORMALISM 

We start our study with some general ideas on the 
scattering mechanism of electrons in a NS junction, as 
illustrated in Fig. [1] As we mentioned above, in a NS 
junction an incident electron from the normal side can 
be reflected by the superconductor as an electron or as a 
hole. The electron and hole wave functions ijj e and ip , re- 
spectively, satisfy the Bogoliubov-De Gennes equation 3 ^ 
given by 



H A 

A* —H. 



(i) 



where Hq is the single electron Hamiltonian and A is 
the potential which couples the electron and hole wave 
functions. We assume that Hq is independent of spin- 
orbit interactions. On the other hand, far from the NS 
boundary the superconducting gap A goes to zero in the 
normal-metal region, whereas in the superconducting re- 
gion the pair potential has an amplitude A and phase <f>. 
Near the NS interface, however, A depends on the details 
of the NS junction; assuming that the width of the NS 
interface is small compared to the superconducting co- 
herence length and neglecting any pairing interaction in 
the normal region, the pair potential is usually modeled 
as a step potential: A(r) = Aoe l ^9(a;), Q(x) being the 
step function. Now we consider that the energy of an in- 
cident electron is in the energy gap (E F < E < E F + A , 
E F being the Fermi energy) where there are no propagat- 
ing wave solutions in the superconducting region, i.e., the 
electron is reflected as a hole or electron. Let us denote 
by r£ e the reflection matrix describing the Andreev re- 
flection. Within the Andreev approximation, where only 
Andreev scattering processes are taking into account and 
with the above assumptions for A, a solution of Eq. (Q]) 
in the normal region can be written a a 18 i 25 



e 
\j) h 



j i exjp[ik^x]4> n (y,z,E F + e) 
+ (®a ) exp[ik^x\(/) n (y,z,E F - e) . (2) 



In this expansion, 4> n (y, z) is the transverse wave func- 
tion that carries a current in the +x direction and n la- 
bels the transverse mode-open channel-, while the wave 
number k% h is given by fcjj = ^/2m/ft 2 V 'E F — E n + e 
and k% = v / 2m/h 2 VE F - E n - e. 

In the linear-response limit, i.e., zero bias voltage (e = 
0), we have fe^j = In this limit, using the continuity of 
the wave function and its derivative at the interface, the 
reflection matrix r^ e is given by r^ e = i cxp (— «</>)1jv and 
r e/i = * ex P where ljv is the N x N unit matrix. 
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FIG. 1: Illustration of a normal-metal-superconductor (NS) 
junction. The coefficients (a e , ah) and (b s ,bh) are related by 
the scattering matrix S A , see Eq. Q. 



A. Andreev scattering matrix 

As we have mentioned, we consider attached leads sup- 
porting one open channel, therefore from now on we spe- 
cialize to this case. Thus, the 2x2 scattering matrix 
for the Andreev reflection which relates the coefficients 
(a ei a h ) and {b e ,b h ), see Fig. [fl is given by 



ah 



with 



where "~ 



1 eh 



(3) 



(4) 



he = iexp(-i(; 



and r£ h 



-(rtr 



B. Schottky barrier 

We also consider the potential barrier created at the 
NS interface due to the mismatch of the conduction bands 
of the two different regions. We model this Schottky bar- 
rier as a planar barrier characterized by its transparency 
r. It is of relevance to consider the effect of this potential 
on the transport since, as we will see later, the conduc- 
tance and therefore its statistical properties are affected 
depending on the strength of T. 

On the other hand, since in the normal region electrons 
and holes are not coupled, the scattering matrix of the 
Schottky barrier S can be written as 



S* = 



S» 

o si 



(5) 



with the scattering matrix for electrons 5f given b y 31 i 32 



(6) 



while the scattering matrix for holes S% satisfies S% 



C. Scattering matrix of the quantum dot 

We now introduce the polar representation for the scat- 
tering matrix S D associated to the quantum dot. For 
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electrons, the 2x2 S^-matrix is written as 



/re' 



re" 



(7) 



N 



where r and the phases 9, 9', ip, and ip' are independent 
variables. In presence of time reversal symmetry, is 
a symmetric matrix and therefore 9 = 9' and ip = As 
before, since in the normal region the electrons and holes 
are not coupled, the complete scattering matrix for the 
quantum dot can be written as 



S 1 







(8) 



with 



III. TWO-TERMINAL CONFIGURATION 

First we consider a two-terminal setup. As shown in 
Fig. [21 a normal lead is attached to a normal quantum dot 
while a second lead with a superconducting region and 
a Schottky barrier of transparency T is also attached. 
Both leads support one open channel. In a similar two- 
terminal setup and for a large number of channels in the 
leads, the average and variance of the conductance have 
been studied in Ref. flol. 



A. Total scattering matrix and conductance 

We are ready to calculate the full scattering matrix 
of the system, Fig. [21 since we have introduced already 
the scattering matrices associated to the components of 
our hybrid structure in the previous section, i.e., we only 
need to combine properly the scattering matrices S A , S , 
and S D from Eqs. (@}, ©, and ©, respectively. 

We combine first the scattering matrices associated to 
the Andreev reflection S A and Schottky barrier S s . The 
resulting scattering matrix S SA is given by 



S ^ = r s +t' s S A (l-r' s S A )~ 1 t s , 



where we have defined the matrices 

rS \ lS ( € 



r 



rf 
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(9) 



(10) 



From Eq. JBJ), rf = rf = VI - F, ff = iVf, and tf 
— iyT 7 , we obtain 
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<->■' e T t e 'eft'/, ' he 1 ' 



+ iS r A +S 
6 /i ' he b e 



f iS r A f S 
b e ' eh b h 



A r S i 4-iS r A r /S r fAj-S 
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(11) 



Schottky barrier — - 



FIG. 2: Two-terminal setup. A normal-metal lead (left) and a 
superconducting lead (down) are attached to a quantum dot. 
Each lead supports one open channel. A Schottky barrier of 
transparency T is introduced in the superconducting lead. 



with A = 1— r' h s r A e r'S rf h . We now combine the above ex- 
pression for S with the scattering matrix of the quan- 
tum dot S D . The resulting scattering matrix S is given 
by 



S = r u +t' u S SM (l-r' D S SB ) ~t L 
After some algebra we write the total matrix S as 

S ■ 

where 

s ee = 2Ce 2 ^ +e '> [(r - 2)VT 



See S e fi 
She Shh 



(12) 



(13) 
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with 



VI - T ((2 - r) cos (tp + i>') + ir sin (-0 + ip')) ' 
_iCr e *(*+*-*'+«-«')lY , 
_iC f e -«( , H-V'-V''+e-9')rT , 
2Ce- 2 W)[(r-2)VT^ : + 
Vl -T((2 - r) cos (V> + ^') - ir sin (0 + V')) 



(14) 



C = [(r - 2)(r - 2) - 4V(r- l)(r- l)cos(^ + 0')]" 1 • 

Finally, the dimcnsionless conductance g due to the 
presence of the NS interface in our two-terminal setup is 
given, in terms of elements of the total scattering matrix, 



G 

9= -77- = 1s he s\ e , 

(jo 



(15) 



where Gq{= 2e 2 /h) is the conductance quantum. 

B. Conductance distribution 

We now calculate the distribution of the conductance 
by assuming that the statistics of the scattering matrix 
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of the quantum dot S D follows one of the circular ensem- 
bles, namely, the circular orthogonal ensemble ((3=1) or 
the circular unitary ensemble ((3 = 2). 

From Eq. (fTS)) , the distribution of the conductance is 
given by 



PW( 9 ) 



9 



s he sl e )d^(S D ), 



(16) 



with invariant measure^ 



d^\S D ) 



dr dO dip 
2^/t 2?r 2tt 



for (3 



d9d^d9'dvb' 

dr for (3 = 2 

2ir 2tt 2tt 2tt 



(17) 



In the following subsections we provide analytical ex- 
pressions for the distribution of the conductance for the 
two-terminal configuration by the evaluation of the inte- 
grals in Eq. (Tl7Jl) . for the cases (3 = 1 and (3 = 2. 




FIG. 3: (Color online) Distribution of the conductance 
P (g) for the two-terminal setup for different values of the 
transparency V. Good agreement is observed between theory 
[Eqs. (|20[) and (|21[) ] and the corresponding numerical simula- 
tions (histograms). 



1. Time-reversal-symmetric case, (3 = 1 

We recall that for the orthogonal symmetry, S D and 
are symmetric scattering matrices and therefore 9 = 
9' and i/j = ip' in the polar representation of Eq. ([7]) . 
Thus, from Eq. flTiJ) we find 



s h e s he 



r t 



(r - 2)(t - 2) - V(r-l)(r-l)cos(2^) 



(18) 

Before considering the general case of arbitrary T let us 
consider the particular case T = 1, i.e., maximum trans- 
parency of the Schottky barrier. In this case Eq. l|T8j) 
reduces to 



s heS he 



(2-rf 



for r = 1 



(19) 



Thus from Eqs. (|T5|l and ([To]) , we obtain the simple ex- 
pression for the distribution 



P [1 \g) = 



i 



2[^g(V2 + ^g)] 



3/2 



for T = 1. 



(20) 



which we have already normalized. 

We now study the case of arbitrary T. Therefore we 
need to perform the integrals in Eq. (fl6|) with S] le s\ e 
given by Eq. (fig)) . The calculation is direct but somewhat 
lengthy. We finally get 



P [1) l9) 



1 



r g 3 / 2 



VB 2 - 4AC - B 
C 3 

[E(x, to) — F(x, to)] 



1/2 



(21) 



where 



A = 16(1 - T) - 4(r - 2) 
B 



C 



4(r - 2y + i6(i - r + 4 A / - (2 - r) r 

F) 2V 2 

2,/-r(r-2)-(r-2) 2 , 

V 9 9 



B + VB 2 - 4AC 
B - ^W^tW ' 

iarcsinh 



2Ct 



B + VB 2 - AAC 



F(x, m) and E(x,m) in Eq. (|2"Tj) are the elliptical inte- 
grals of first and second kind, respectively, which have to 
be evaluated at the limits r + and r_ given by 




(22) 

We have verified our predictions for P^^g), Eqs. (|2"0|) 
and (|21[) . by comparing with numerical simulations of a 
quantum-chaotic dot (in all our numerical simulations we 
have considered a Bunimovich stadium geometry whose 
classical dynamics is chaotic). In Fig. [31 we show both 
numerical and analytical results for T = 0.2, 0.5, and 
1. We observe good agreement between theory and the 
corresponding numerical simulations. 



2. Broken time-reversal symmetry, (3 = 2 

We now consider the circular unitary ensemble. We 
recall that this symmetry class is relevant for situations 





FIG. 4: (Color online) Distribution of the conductance 
p( 2 \g) for the two-terminal setup for different values of the 
transparency T. Good agreement is observed between theory 
[Eqs. (|24[) and (|25fl] and the corresponding numerical simula- 
tions (histograms). 



FIG. 5: Three-terminal setup. Two normal-metal leads (1 and 
2) and a superconducting lead (3) are attached to a quantum 
dot. Each lead supports one open channel. A Schottky barrier 
of transparency T is introduced in the superconducting lead. 

become statistically favored as the Schottky barrier be- 
comes transparent. We anticipate that the effect of the 
Schottky barrier on the conductance is somewhat differ- 
ent in the three-terminal configuration that we study in 
the following Section. 



where time-reversal symmetry is absent. From Eq. (|14[) 
we have 



S he s he 



2T 2 r 



2^2 



(r - 2)(r - 2) - V(r-l)(r-l)cos(V> + V) 

(23) 

As in the previous subsection we start with the special 
case r = 1. Thus, from Eqs. ((HJ), (pH, and ® the 
distribution of conductance is given by 



P (2> (9) 



1 



for T = 1. 



(24) 



The general case of arbitrary T involves a direct but 
lengthy calculation, as in the (3 = 1 case. For (3 = 2, 
however, the distribution gets a simpler expression. Af- 
ter some algebra we find 



P {2) (9) 



[r(V2^2) + 4] r 2 
Vg[V2r + V5(2-r)] : 



(25) 



In Fig. H we plot P (2) (ff), as given by Eqs. ([24]) and ([25]) . 
and compare with numerical simulations of a quantum- 
chaotic dot. We choose T = 0.2, 0.5, and 1. As for the 
case (3 = 1, for broken time-reversal symmetry a good 
agreement between theory and numerics is also observed. 

To conclude this Section we notice that, for the two- 
terminal configuration studied here and for both symme- 
tries {(3 = 1 and (3 = 2), the distribution of the conduc- 
tance becomes broader as the value of the transparency 
of the Schottky barrier is increased, see Figs. [3] and |U 
This broadening of the distribution implies that large 
values of g (values of g at the tail of the distribution) 



IV. THREE-TERMINAL CONFIGURATION 

We now consider a more complex situation. As shown 
in Fig. El two normal leads are attached to a normal quan- 
tum dot while a third lead with a superconducting mate- 
rial and a Schottky barrier is also attached. All the leads 
support one open channel. As in the previous Section, we 
calculate the total scattering matrix of the hybrid struc- 
ture. The present case is, however, algebraically more 
involved than the two-terminal configuration. 



A. Total scattering matrix and conductance 

We define a scattering matrix S that relates the incom- 
ing and outgoing wave amplitudes in the normal leads 
(leads 1 and 2 in Fig. [5|): 



fb \ 




bl 


- s [i 


U 2 J 





(26) 



where S has the following structure 



S = 



( 



'ii 

„ee ee 
'21 s 22 



'11 



'12 a ll 
e eh 
s 21 
hh 

'12 *11 



eh \ 



-, 2 2 
.hh 
'12 



„he „he „hh ~hh 



(27) 



Here, the subscripts denote the leads. See Fig. [5] We 
recall that we are working at energies where there are 
no propagating waves in the superconducting region. In 
order to obtain the explicit expressions for the elements 
of the scattering matrix S, we need to combine the scat- 
tering matrix of the superconducting region [Eq. (j4|] , the 
Schottky barrier [Eqs. © and ©], and the quantum dot. 
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The scattering matrix associated to the quantum dot S D 
can be written in two blocks, as in Section III Ci one for 
the electrons and the other for the holes. Thus, for the 
three-terminal setup, we write for electrons 



t 



31 °32 



'12 l l 

'22 '2 
t 



(28) 



' 33 



with a similar expression for the holes. Once we have 
introduced the scattering matrices associated to the dif- 
ferent parts of our setup, we can follow the same steps of 
the previous Section to obtain the total scattering matrix. 
On the other hand, however, the dimensionless conduc- 
tance g for our three-terminal setup is given by^ 



respectively, arc simplified to 



'23' 33 '31 



'21 



'21 



' 33' 33 



t h t 



23 '31 



e h 
'33' 33 



Therefore the conductance is given by 



j-e „h j-e 
t 23 , 33 t 31 


to 


jh j-e 
'23'31 


1 + r 33 r 33 




1 + r 33 r 33 



2. Case T = 



(32) 
(33) 

(34) 



„ee |2 



-eh 1 2 



(29) 



Thus, we need only the elements s|f and s|j of the total 



scattering matrix S, Eq. (|27|) . in order to calculate the 
conductance. After some algebra we find for these two 
scattering-matrix elements 



For the case of total reflection at the Schottky barrier 
the conductance is given by 



*21 + 



j-e j-e 
t 23' , 31 



1 + r 



33 



(35) 



where 



1 

A 



A% + -C e 2 r^ h C h r±A e 



1 

A 



■ °2 r he A 



(30) 
(31) 



Clearly, there is no electron-hole scattering processes in 
this case. 



B. Conductance distribution 



A 

A', 



1 ~ r he^ er ebp h 

j-e r iSj.e 
'23' e s 31 



*21 + 



1 - ryr'S 



CS = t 



23 



r lS r e j-S 
j_S _|_ r e r 33'e 



1 - r 33 r e S 



A' 



,lS,e _ 
l e '31 



' 33 ' h 

j-iS„e -iS 
l e '33 c e 



1 - n3 r e 



IS 



C h and C! 2 l in Eqs. and (5TJ) are defined as C e and 
Cf, respectively, by replacing e — » /i. We can observe 
that if the third lead were absent, the expression for con- 
ductance, Eq. (|29|) . is reduced to the usual Landaucr- 
Biittiker formula g = \t\^ = |*2i| 2 - 

For the particular values of T = and 1, the expres- 
sions for the elements s|f and s|i are simplified. There- 
fore for these two special values of T, it is useful to show 
the expressions of s|i an d s|i si nce they exhibit explicitly 
the different scattering processes among the three leads 
attached to the quantum dot that contribute to the con- 
ductance. 



Case r = 1 



For the case of zero reflection at the Schottky barrier 
the expressions for s 2 \ and in Eqs. (|3"U|) and (|3"Tj) . 



In order to calculate the distribution of conductances 
we need to average over the orthogonal or unitary ensem- 
bles ((3=1 and (3 = 2), as in the previous Section. Sec 
Eq. (|16j) . In the three-terminal configuration, however, 
the expression for the conductance [Eqs. ([2"9")) - (prr]) ] is 
more complicated, even at the values of the transparency 
r = and 1, where the expressions for the scattering- 
matrix elements sff and s|i are simpler. Therefore we 
are not able to provide analytical expressions for P(g), 
as for the two-terminal setup. Thus we proceed with 
our analysis by generating numerically an ensemble of 
conductances from Eqs. (f2"9")) - (|3"2"l) . where the elements of 
the scattering matrix associated to the quantum dot are 
generated according to the orthogonal and unitary en- 
sembles. 

In presence of time- reversal symmetry ((3 = 1), the 
distribution of the conductance at different values of the 
transparency is shown in Fig. [SI At zero transparency, 
the distribution is described b y 36 ' 37 P^\g) = \j\f%g 
(solid line in Fig. [6]), which is the known result for 
a quantum-chaotic dot with two normal single-channel 
leads, i.e., there is no effect of the superconducting region 
on the conductance, as one might expect for a nontrans- 
parent Schottky barrier. For T > 0, however, P^\g) is 
affected. From Fig.[6]we can see that for intermediate val- 
ues of the transparency (r = 0.2 and T = 0.5) the distri- 
bution probability is finite for g > 1 , whereas for null and 
complete transparency (r = and T = 1, respectively) 
the probability is confined in the interval < g < 1. 
Therefore multiple scattering in the third lead, where 
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FIG. 6: Distribution of the conductance P^(g) for the three- 
terminal setup for different values of the transparency V. No- 
tice that P {1) (g > 1) > for < T < 1, see right panels. The 
solid curve for T = is given by P (<?) = l/y/4g. 
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FIG. 7: Distribution of the conductance P^ 2 \g) for the three- 
terminal setup for different values of the transparency V. As 
for (3 = 1, here P (2) (.9 > 1) > for < T < 1, see right panels. 
The horizontal line for F = corresponds to P^ 2 '(g) = 1. 



the Schottky barrier is present, produces larger values of 
conductance. 

When time reversal symmetry is broken ((3 = 2) the 
distribution of the conductance is quite different from the 
(3 = 1 case. In Fig. [7|we show P^ 2 \g) at different values 
of the transparency. First we notice that for r = 0, 
P^ 2 \g) reduces to the known resul t 36 ' 37 P^ 2 \g) = 1 
(shown as the solid line in Fig. [7]). At intermediate val- 
ues of the transparency of the Schottky barrier (r = 0.2 
and r = 0.5), we observe again that multiple scattering 
in the superconducting lead can produce larger values of 
conductance than for the cases T = and T = 1, i.e., the 
distribution P' 2 ' (g) develops a tail for g > 1 at interme- 
diate values of T. 



V. SUMMARY AND CONCLUSIONS 

We have studied the effects of a superconducting lead 
attached to a normal quantum-chaotic dot on the sta- 
tistical properties of the conductance using random ma- 
trix theory and a scattering approach to the transport in 
normal-metal-supcrconducting (NS) hybrid structures. 
We have calculated the complete distribution of the con- 
ductance in the presence and absence of time-reversal 
symmetry (/3 = 1 and (3 = 2, respectively) and for two 
different configurations: two and three-terminal setups. 
The conductance distribution has been analyzed as a 
function of the transparency T of the Schottky barrier 
which is modeled as a planar barrier formed at the NS 
interface. The importance of considering the effect of a 
potential barrier formed at NS junctions has been already 
recognized in transport experiments^ however, opposite 
to the one open channel case that we study here, in these 
experiments a large number of channels are open. 

For the two-terminal configuration we provided analyt- 
ical expressions for the distribution of the conductance. 
We found that the distribution becomes more and more 
concentrated at g — as the value of the transparency 
is reduced, for both symmetries (3 = 1 and (3 = 2. This 
behavior might be expected considering that the conduc- 
tance is due to the scattering at the NS interface (An- 
dreev reflection), i.e., at zero transparency of the barrier 
the situation is reduced to a normal quantum dot with a 
normal lead attached where only reflection processes at 
the lead are present. 

For the three-terminal configuration, the analytical ex- 
pressions for the total scattering matrix and the conduc- 
tance are more complicated than for the two-terminal 
case. We thus obtain the distribution of the conductance 
from an ensemble of conductances generated numerically 
according to our analytical expressions for the elements 
of the total scattering matrix. We found that in presence 
and absence of time-reversal symmetry, conductance val- 
ues larger than unity are possible at intermediate values 
of the transparency of the Schottky barrier. As a con- 
sequence, the distribution of the conductance is broader 
than for the cases of zero transparency (r = 0) and an 
ideal NS interface (r = 1), i.e., multiple scattering in 
the superconducting region with the Schottky barrier in- 
creases the possibility of larger values of the conductance. 

To conclude, we have seen that the presence of a Schot- 
tky barrier in the lead, which is formed in real experi- 
ments due to the mismatch of the conduction bands of 
the normal-metal and the superconductor, has a signifi- 
cant effect at the level of the distribution of the conduc- 
tance. The progress in the fabrication of nanoscale het- 
erostructures has made possible the experimental study 
of statistical properties of the conductance through such 
devices^ These experiments have pointed out the rel- 
evance of non-perfect contacts on the study of conduc- 
tance. We thus think that the effects on the statistics of 
the conductance through quantum dots due to the pres- 
ence of Andreev scattering and potential barriers at the 
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leads, as those described here, are of experimental inter- 
est. 
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